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Abstract 

Within the complete heavy quark effective field theory (HQEFT), the QCD 
sum rule approach is used to evaluate the decay constants including I/tuq 
corrections and the Isgur-Wise function and other additional important wave 
functions concerned at I/tuq for the heavy-light mesons. The number of 
unknown wave functions or form factors in HQEFT is shown to be much 
less than the one in the usual heavy quark effective theory (HQET). The 
values of wave functions at zero recoil are found to be consistent with the 
ones extracted from the interesting relations (which are resulted from the 
HQEFT) between the hadron masses and wave functions at zero recoil. The 
results for the decay constants are consistent with the ones from full QCD 
sum rule and Lattice calculations. The l/mq corrections to the scaling law 
/m ~ F/y/rriM are found to be small in HQEFT, which demonstrates again 
the validity of 1/mg expansion in HQEFT. It is also shown that the residual 
momentum v ■ k of heavy quark within heavy-light hadrons does be around the 
binding energy A of the heavy hadrons, which turns out to be in agreement 
with the expected one in the HQEFT. Therefore such a calculation provides a 
consistent check on the HQEFT and shows that the HQEFT is more reliable 
than the usual HQET for describing a slightly off-mass shell heavy quark 
within hadron as the usual HQET seems to lead to the breakdown of 1/mQ 
expansion in evaluating the meson decay constants. It is emphasized that the 
introduction of the 'dressed heavy quark' mass is useful for the heavy-light 
mesons (Qq) with thq >> A >> m q , while for heavy-heavy bound states 
(ipifa) with masses mi, mi >> A, like bottom-charm hadrons or similarly 
for muonium in QED, one needs to treat both particles as heavy effective 
particles via \jm\ and 1/mi expansions and redefine the effective bound 
states and modified 'dressed heavy quark' masses within the HQEFT. 
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I. INTRODUCTION 



It has been seen that the heavy quark effective field theory (HQEFT) with keeping both 
quark and antiquark fields Jl| can provide a consistent description on both exclusive [0 and 
inclusive @,|] decays of heavy hadrons. The extracted values of |V^,| from both exclusive 
and inclusive decays have shown a good agreement. The lifetime differences among bottom 
mesons and hadrons can also be well understood within the new framework of HQEFT. 
Especially, it has been noticed that at zero recoil, I/tuq corrections in both exclusive and 
inclusive decays are automatically absent without imposing the equation of motion iv-DQ v = 
when the physical observables are presented in terms of heavy hadron masses, this is the 
main point that differs from the usual heavy quark effective theory (HQET) or the usual 
heavy quark expansion. In the usual framework, 1/m.Q corrections in the inclusive decays 
are absent only when the inclusive decay rate is presented in terms of heavy quark mass 
{ m o) @ rather than the heavy hadron mass (ra#), the situation seems to be conflict with 
the case in the exclusive decays where the normalization is given in term of heavy hadron 
mass || . Such an inconsistency in the usual HQET may be the main reason that leads to 
the difficulty for understanding the lifetime differences among the bottom hadrons. These 
observations indicate that the contributions and effects from antiquark fields should play a 
significant role for understanding hadronic structures and can become important for certain 
physical observables. 

Our basic point of considerations is based on the physical picture that a heavy hadron 
(Qq) containing a single heavy quark (Q) and satisfying the condition among the heavy 
quark mass itlq, the light quark mass m q and the binding energy A 

rriQ » A >> m q 

may be regarded as a 'dressed heavy quark' with an off-mass shell by an amount of binding 
energy A. In this case, the usual quark-hadron duality should be extended to a 'dressed- 
heavy quark' - 'heavy-light hadron' duality. Thus a more reliable heavy quark expansion 
for heavy-light hadron systems should be carried out in terms of the "dressed heavy quark" 
mass defined as 

tjiq = m Q + A = m H - 0(l/m Q ) = Jmv^m H 

with rriH the heavy hadron mass. Before proceeding, we would like to point out that this 
picture cannot naively be applied to the hadrons with 

m Q ,m q » A 

This is because such hadrons must be treated as heavy-heavy hadrons, like B c meson system. 
Similarly, the muonium (fie) must also be treated as heavy-heavy bound state system in 
QED. In general for heavy-heavy system (^1^2) with masses m 1 ,m 2 >> A, one should 
make expansion for both heavy particles in the bound state and redefine the effective fields 
and effective bound states. 

It has been shown |^@] that the new framework of HQEFT with keeping the antiquark 
fields enables us to describe such an off-mass shell heavy quark within hadrons as the 1/mg 
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corrections to the transition matrix elements automatically vanish at zero recoil. Of partic- 
ular, the number of transition form factors invloved in the HQEFT is much less than the one 
in the usual HQET. In addition, the new framework of HQEFT also results in some interest- 
ing relations among meson masses and transition formfactors at zero recoil, which enables 
us to extract the important transition form factors from the known heavy meson masses. 
It is those special features that allow us to check self-consistently the validity of the new 
framework of HQEFT when applying it to the heavy-light hadron systems. This can simply 
be carried out by comparing the values of transition formfactors and residual momentum of 
heavy quark in hadron, which are extracted from the relations between the meson masses 
and formfactors, with those obtained from directly evaluating the corresponding hadronic 
matrix elements by using some reliable approaches. In this paper, we are going to adopt 
QCD sum rule approach for a practical calculation. Our paper is organized as follows. In 
Sec.[H], we briefly review the heavy quark expansion and present a general description on 
form factors concerned up to order I/tuq, their corresponding hadronic matrix elements in 
heavy meson weak decays and weak transitions between two heavy mesons are expilcitly de- 
fined. In Sec.[TTI|, we apply the QCD sum rules to two point Green functions of heavy quark 
currents within the framework of HQEFT, and calculate three form factors F, G\ and G2 
as well as two composite factors gi and concerned in the decay and coupling constants 
of heavy pseudoscalar and vector mesons. It also enables us to extract a reasonable value 
of the binding energy A. In Sec.[IV], we investigate three point Green functions of heavy 
quark operators and apply the QCD sum rule to evaluate the Isgur-Wise function and some 
additional wave functions appearing at the order of I/tuq. It is also shown that the residual 
momentum of the heavy quark within the hadron is truely around the binding energy, i.e., 
iv ■ D ~ v ■ k « A is found to be a good approximation in simplifying the evaluation of the 
hadronic matrix elements. In Sec.[V|, we add the two-loop perturbative contrinutions to the 
two-point collerator for the purpose of seeing the importance of QCD radiative corrections. 
A brief summary with some remarks is presented in Sec.[VT]. 



II. HADRONIC MATRIX ELEMENTS IN HQEFT 



It has been shown that by decomposing the original heavy quark field in QCD into 
effective quark field and antiquark field [0] and integrating out both the antiquark field 
and the small components of quark field, we arrive at an effective theory for only the large 
components of quark field Q+, its effective Lagrangian has the following form |||3| 

r(++) _ r(0) , r(l/m Q ) 
^eff ~ ^eff + L eff 

= Q+i%Qt + -^-QtW±) 2 Q+ + 0{\), (2.1) 

m Q m Q 

where is the heavy quark effective field in this new framework of HQEFT, its momentum 
is k = Pq — rrigv with Pq the momentum of the original heavy quark field in QCD. The 
operators Ipu and Ip± are defined as 

f g(x) D» <p(x) = - f g(x)Di><p(x), 
l % =i>{v-D), (2.2) 



3 



with v M an arbitrary four- vector satisfying v 2 = 1. In general, a mass dimension parameter 
A may appear in the Lagrangian depending on the redefinition of heavy quark effective field 
0, here we take A = for the convenience in comparison with the convention in the usual 
heavy quark effective theory(HQET). 

The decay constants of a heavy pseudoscalar meson P and a heavy vector meson V are 
defined by 

< Q\q~J»j5Q\P(v) >= ifpmpVp, 

<0\q^Q\V(e,v) >=f v m v e lx (2.3) 

with Q the original quark field in QCD. Here e M is the polarization vector of a vector meson. 
\P(v) > and \V(e, v) > are the pseudoscalar and vector meson states in QCD. 

In the new framework of HQEFT, current operators composed by one heavy quark and 
one light (or heavy) quark can be expanded in terms of operators in the effective theory 

qTQ -> qTQt + J-gT^L(i]p ± ^Q+ + 0(l/m 2 Q ), 
2m Q iip\\ 

qtq -> g;trg+ + J-Q>+r^-(iij) x ) 2 Qt 

2m Q ilp\\ 

1 Q:tH^) 2 ^rg+ + o(iA4). (2.4) 
-i % 



2m Q > 



Correspondingly, it is useful to introduce an effective heavy hadron state \H V > for 
exhibiting a manifest heavy quark spin-flavor symmetry 0. Such an effective hadron state 
in HQEFT is related to the hadron state \H > in QCD via P|| 

< H'\Q'TQ\H >= _ < H^J^e'J^ff \H V > . (2.5) 



V m H'm H " " y/A HI A H 

with 

A H C) = m H(') ~ m Q(')- ( 2 -6) 
The renormalization condition for \H V > is given by 

< H V \Q+YQ+\H V >= 2Av» (2.7) 

with 

A = A H - 0(l/m Q ) = lim A^ 

being a heavy flavor independent binding energy that reflects the effects of the light degrees 
of freedom in the heavy hadron. Obviously, the above normalization condition preserves 
spin-flavor symmetry. We would like to address again that above renormalization is only 
applicable for the heavy-light hadrons with uiq » A » m q . For heavy- heavy bound state 
system, i.e. mi, m 2 » A, such as bottom-charm meson B c and muonium (fie), one needs to 
redefine the effective bound states by considering 1 / uiq expansion for both heavy particles. 
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In the heavy quark expansion, 1 / mg order corrections to the hadronic matrix elements 
arise from not only the current expansion Q2.4Q but also the effective Lagrangian, i.e., inser- 



tion of L^jj 1 ^ into the matrix elements. By including all these corrections, we arrive at the 



following results for the hadronic matrix elements of the currents in ( F2.4|) 

< 0\qTQ\M >^< 0\qTQ+\M v > --L < 0\qT±-{ilp ± )*Q+\M v > +0(l/mj 



X < M'\Q'YQ\M >^< M^,\Q' v tTQ+\M v > < M^lT— (iIp ± ) 2 Qt\M v > 

\ m M >m M 2m Q ilp\\ 

1 < M^\Q'+(-i $ ± ) 2 -^rQt\M v > +0(l/m|). (2.8) 



2m Q/ 

The relevant hadronic matrix elements may be parameterized as follows 
<0\qTQt\M v >=^Tr[TM], 

< 0\qT—^P + (D ± ) 2 Q+\M v >= -FGiTrFM], 

< 0\qT—^P + ^a al3 F^Q+\M v >= 2FG 2 Tr [ia a pTP + ^ a al3 M] = -2FG 2 d M Tr[TM], 

< M' V ,\Q'+TQ+\M V >= -C(y)Tr[M'TM], 

< M^\Q'+V^-P + DiQt\M v >= —Ki(y) ^Tr[M'TM), 

iv ■ D A 

< M'JQ'+T—^P + l -a aP F^Q+\M v >= iTr[K aP (v,v')M'TP + l -a^M] } (2.9) 



where y = v ■ v '. Ai(v) is the spin wave function 

M( v ) — v^AP I — ^ or P seu( ioscalar meson P 
+ 1 for vector meson V 



(2.10) 



and 



^ _ I 3, for pseudoscalar meson P (2 11) 

1—1, for vector meson V, 

where P, G\ and G2 ar e constants, and £(y) and /«i(y) the Lorentz scalar functions. Ac- 
tually, £(?/) is the well-known Isgur-Wise functions. The Lorentz tensor K a p(v,v') can be 
decomposed into 

Kap(v,v') = iK 2 {y)a a p + K^vWalp - «j7o) ( 2 - 12 ) 

with ^2(2/) and Ks(y) being the Lorentz scalar functions. 
Combining (|2.3|), ( |2.7|) and fl2.9|) we have 



A 1 



\A^ F ( 1 + ^ (G ' + M »' G2 »^ (213) 
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Thus the ratio between the vector and pseudoscalar meson constants is given by 



fprrij!, Ay rriQ 

As is known, the normalization of the Isgur-Wise function at zero recoil point is given 
by £(1) = 1. The additional wave functions Ki(y), ^{y) and K^{y) characterize the next-to- 
leading order symmetry-breaking corrections to £. From (|2.9|) and fl2.12|) , it is easily seen 



that the hadronic matrix element at zero recoil is irrelevant to ^3(1). While fti(l) and ^2(1) 
are found to be related to the meson masses 0|| 

Am = fn M — rriQ = A — (— __)(k 1 (i) + d M K 2 (l)). (2.15) 

m Q Zm Q 

It is easy to check that the leading order contributions to the hadronic matrix elements in 
the HQEFT, which are characterized by the decay constant F and the Isgur-Wise function 
£ defined in ( |2.9| ), are the same as the ones in the usual HQET. Nevertheless, to the next- 
to-leading order, differences occur between the two frameworks. In the usual HQET, the 
transition matrix elements between two heavy mesons are parameterized by six functions 
denoted as £j and %% (* = 1> 2, 3). Here £1, £2 and £3 arise from the current expansion, and \u 
Xi and X3 from the insertion of the effective Lagrangian into the hadronic matrix elements. 
All these six quantities are functions of the recoil parameter y — v • v'. In addition, there are 
two more parameters Ai and A2 which appear in the l/mq order corrections to heavy meson 
masses. Unlikely, in the new framework of HQEFT, it has been noticed that the evaluation 
of the hadronic matrix elements is greatly simplified when the antiquark contributions are 
included As a consequence, all transition matrix elements concerned at 1/ttlq order 

can be characterized by only three wave functions K\{y), ^(y) and K%(y). Furthermore, the 
1 / mq order corrections to meson masses were found to be naturally related to their values at 
zero recoil, i.e., and ^(l)- In other words, instead of the six wave functions £i(y) and 

Xi(y) (i— 1)2,3) and two parameters Ai and A2 in the usual HQET, we only need to evaluate 
three wave functions Ki (i = 1, 2, 3) in the new framework of HQEFT up to the order I/tuq. 
At zero recoil, only two parameters fti(l) and ^2(1) are relevant. Similar comments hold for 
the corrections to meson decay constants. In the HQEFT one encounters only two constants 
G\ and G2 at 1/tuq order. 

Within the new framework of HQEFT, two important parameters /«i(l) and /^(l) have 
been extracted from meson mass spectrum in refs. Jl . To have an independent check for the 
HQEFT, it would be very useful to evaluate these two wave functions directly by a field- 
theoretical method within the framework of HQEFT. In the following sections, we shall 
present an QCD sum rule study for these form factors. 



III. QCD SUM RULE CALCULATION OF F, G x AND G 2 

QCD sum rule approach has widely been used to calculate hadronic matrix elements in 
QCD and has also been applied to effective theories of QCD. It turns out to be a powerful 
analytic approach to estimate non-perturbative effects. The basic idea of QCD sum rule 
formalism is to study the analytic properties of correlation functions, and to treat the bound 
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state problems in QCD from quark-hadron duality considerations. So that one could start at 
short distance physics and moves to large distance physics where confinement effects become 
important and resonances emerge as a reflection of confinement. Here we shall briefly outline 
the main steps of sum rule treatments. In general, the fixed point gauge for the background 
fields is used in calculating Feynman diagrams. 

In order to evaluate the parameters F, G\ and G2 involved in the meson decay constants, 
we consider the following two-point correlator 

U(u) = 1 J d*xe iPx < 0\(qT MI Q) {x) , (Qr M g) (0) |0 >, (3.1) 



where Tm has appropriate Lorentz structure so that the two currents in ( |3.1| ) interpolate 
the heavy meson of interest. It is convenient to choose || 



M 



—27 s , pseudoscalar meson P , . 

7m — v tn vector meson V 



The total external momentum in (3J.) is P = rriQV + k with the momentum k in (3JJ being 
the residual momentum of the heavy quark. The correlator II is an analytic function of 
2v ■ k + k 2 /rriQ with discontinuities for its positive values. Here fc M = fcj. + k£ with fcj. and 
k^ = (v- k)v^ being the transverse and longitudinal part of k separately. Particularly, under 

k 2 

the definition 00 = 2v ■ k H — — , one has 

m Q ' 

1 2 2 

2vk+ = lu + ^— + Oil/ml). (3.3) 

In eq. ( p.l| ) we have represented n as an analytic function of the variable u. 

Phenomenologically, the two-point function H(uj) can be written as the sum of three 
parts: a pole contribution from the ground state mesons associated with the heavy-light 
currents; a dispersion integral over a physical spectral function; and subtraction terms, 
namely, 

tt < \ ^<0\qT M Q\M><M\QT M q\0> 

UphenM - -(2J P2_ m 2 T jf 

pole ^ m M + le 

+ f°° du Pphys ^\ + subtractions, (3.4) 
Juj c v — uj — ie 

where the two matrix elements should be expanded by using eq.( [2.8| ). (J2 po ie) nieans sum- 
mation over polarization for vector mesons. 

On the other hand, using the Feynman rules of the HQEFT, the correlation function in 
HQEFT is evaluated perturbatively in the deep Euclidean region (u A) via 

n(w) = i J d 4 xe ik - x {< 0\(qT M ,Q+) {x) , (Q+r M g) (o) |0 > 

— <0\(qT M ,l 
2m Q iip\\ 



+ 7^ < 0\(qT M ,Q+) {xh (Q+H jPi-f — r A/ g) (0) |0 >} 

2m Q -1 %\ 

+ 0{l/m%) (3.5) 



7 



In lower Euclidean region, the non-perturbative contributions become important. In QCD 
sum rule analysis, the two-point correlator receives contributions not only from the pure 
perturbative ones but also the ones from condensates, which characterizes non-vanishing 
vacuum expectation values of local operators in the operator product expansion(OPE). It 
is thought that by including these condensates, the QCD confinement effects may be ac- 
counted for at the transition from perturbative region to non-perturbative region. Writing 
the perturbative contributions as the form of an integral over a theoretic spectral function, 
the theoretical result for the two-point correlator ( |3.1|) becomes 

Htheoriw) = ( dv Pp ert \ _[_ Y[ NP -f subtractions. (3.6) 
J (is — uj — it) 

A basic assumption in QCD sum rule is the quark-hadron duality. Due to this duality 
one can model the contributions of higher resonance states by the perturbative continuum 
starting at a threshold energy uj c . In other words, we assume p p h ys = Ppert- Equating the 
phenomenological side and the theoretical side, up to the order of 1/tuq one arrives at 

2Tr[T M P + T M ]^[l + + 2d m G 2 )}^^-(l - *%*) 

= Io c dv ^lZ-l) + n ^ + subtractions (3.7) 

with ujm = 2A M . In deriving ( |3.7|) we have used ( |2.9|) as well as the relations 

Tr[f M M(v)]Tr[M(v)r M ] = -2ATr[T M P + T M ] (3.8) 



and 



P 2 -m 2 M u + u) M , 2 



[UJ 



m Q " -M)[l + ^ + 0(l/mi)]. (3.9) 

The relevant Feynman diagrams are plotted in Fig.l. Fig. 1(a) is the lowest order per- 
turbative diagram. For the non-perturbative effects, it is sufficient in the present case to 
consider only the contributions of the quark condensate, the gluon condensate and the mixed 
quark-gluon condensate, which have values (a s = gl/An) 

<qq>^ -(230) MeV 3 ; 

% < qa a pF al3 q >^ -0.8 < qq >; 

a s < FF >= a s < F^Ff >^ 0.04 GeV 4 . (3.10) 

We only keep terms up to order a s for the condensates. Neglecting the light quark mass, 
we obtain the following results at the renormalization scale p ~ 2A 

F 2 [l + — (Gx + 2d M G 2 ) - J_ (^+^g) + A][J_ + ?__(! _ = 

m Q m Q 2A m Q 4m Q uj m - uj - le 2m Q 



, . 1 1 , 3z/ 3 . A 1 or, 4 9-4dj 
dv[- + -](z/ 2 - )+ <qq>[-~ ~ A + -(tt- 



Q 



8tt 2 Jo ArriQ v — lj — ie Avtlq uj Avtlq ty 3oj 9m. 

_ & 1 1 3 d M a 2 163 + 147d M 

+i < qa af} F ap q > — -(- - — + - — + — — — ) 

2u 6 vtlqlo 1 8 12 n oj 6 blbrriQU) 2 

Ot r 1 15 — Ad Mi 

< FF > — + subtractions. (3.11) 
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In order to improve the convergence and suppress the importances of higher-resonance 
states, we apply the Borel operator 

Bp=T 1 ^_„_ > „ W) (-- r (3.12) 

to both sides of (|3.11|) with T = — being held fixed. In the dispersion integral, this 
Borel transformation yields an exponential damping factor which effectively suppresses high- 
resonance contributions. On the non-perturbative terms, Borel transformation enhances the 
importance of low dimension condensates. Furthermore, the subtraction terms in ( |3.11| ) may 
also be eliminated by this Borel transformation. 
The resulting Borel transformed sum rule reads 

F 2 [l + —{Gi + 2d M G 2 ) - -Ltf + d M J$) + A _ ^L ]c -m/t = 
rriQ 2m Q A m Q 2m Q 



8tt 2 Jo ^rriQ 

_ < > (1 + _ , < ^ q > [_L + _1_ ( 3 _ 

+ a -(± - 16 ^ U f' )}- a -<FF>(±- + 1 -^iL). (3.13) 
7i 1 z htbrriQl tt 241 9brriQ 

Besides the I/tuq corrections shown explicitly in ( |3.13| ), the pole energy ujm and the 
threshold energy uo c also receive 1 / mg corrections. We may write them as 

u>m = 2 Am = &m H (wm + c?m^m ) , 

m Q 

uj c = u H (uji + d M u 2 ). (3-14) 

m Q 

With these formulae, eqs. (|2.13|) and (|2.14|) can be rewritten as 

fu = -^={1 + —(9i + 2d M 92)} (3.15) 

and 

fv m V i 8 /o 1K N 

m = 1 g 2 3.16 

fprnj, 2 m Q 

with g\ and g 2 being two composite factors defined as 



gi = Gi 



UJ 



(i) 

M 



4A ' 



(2) 



92 = G 2 -^L. (3.17) 

Now we first consider the sum rule ( |3.13| ) at leading order. For the leading terms, the 
results are the same as the ones in HQET, because there is no difference between the HQEFT 
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and the usual HQET in the limit mg — > oo. We shall not repeat the analysis for the leading 
order sum rule analysis, our numerical results are presented in Fig. 4, where we have used 
a s (2A) ~ 0.34. From the stability of the curves, we are led to the following solutions for the 
parameters 

co = 1.8 ± 0.3 GeV, 
(o) 

-M. = A = 0.53 ± 0.08 GeV, 

F = 0.30 ± 0.06 GeV 3/2 . (3.18) 



We now proceed to the next-to- leading order sum rule analysis in ( p. 13 ). Putting (p. 14 ) 



into ( [3.13j ), and expanding ( |3.13|) in 1/rriQ, we obtain, at 1/tuq order, two sum rule formulae 



which are relevant to the spin-symmetry conserving and violating corrections, respectively. 
These two sorts of corrections are easy to be distinguished because the latter is proportional 
to (1m- To find out the solutions, it is useful to first evaluate the quantities LO\ and uj 2 by 
requiring optimal stability of uj^ and with respect to the Borel parameter T in the 



allowed sum rule windows. Using the central values of lu 0i A and F in (|3.18|) , we then obtain 



for uji, Gi and g± the numerical results plotted in Figs. 5- 10. When the Borel parameter 
T takes the reliable values T = 1 ± 0.2 GeV, one can read off the following solutions 



2 



wi = 1.5 ±0.2 GeV 
u$ = 0.86 ±0.10 GeV 2 , 

d = 0.95 ±0.15 GeV, 

g x = 0.54 ±0.12 GeV, 

to 2 = -0.15 ±0.05 GeV 2 , 
J§ = -0.16 ±0.03 GeV 2 , 

G 2 = -0.09 ± 0.03 GeV, 

g 2 = -0.06 ± 0.02 GeV. (3.19) 

In the usual HQET || two parameters Gi and G 2 were defined for the 1/mq corrections 
to decay constants. The sum rule calculation in || yielded an unexpectedly large value 
for G\ ~ — 4A pa — 2.0GeV, which leads to the breakdown of I/tuq expansion for 

decay constants. In the new framework of HQEFT, however, as can be seen from eqs.( p,15"|) 
and ( |3.16| ), 1/rriQ corrections to the physical decay constants $m and the ratio are actually 
characterized by the composite factors g^. Though G\ in eq.( |3.19| ) is large, — remains small 
enough so that the 1 / rriQ expansion in the new framework of HQEFT appears to be more 
reliable. 

When taking the typical values for the quark masses m;, = 4.8 ± O.lOGeV and m c = 
1.35±0.10GeV, we obtain from ( 3.18|) and (|3.19|) the following values for the decay constants 



of bottom and charm meson without including QCD corrections caused by the running 
energy scale from \x ~ to fi ~ 2A 

/b(2A) = 0.135 ±0.035 GeV, /b*(2A) = 0.147 ±0.034 GeV, 

/d(2A) = 0.246 ± 0.097 GeV, f D . (2A) = 0.308 ± 0.091 GeV. (3.20) 
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Note that all the quantities , cjq, -F\ Wm, k>», Gj (i = 1,2) and $ obtained by QCD 
sum rules are scale dependent. And the results in ( |3.18| ) and (PTT3) are corresponding to 
the values with the renormalization scale at fi ~ 2A m lGeV. So do the decay constants in 

From the results given in ( |3.18| ) and ( |3.19| ), we arrive at the following relations from 
eq. (|3~T6D 



j. 1/2 
, 1/2 

fDm)i 2 



1.10 ±0.03, 



1.31 ±0.07. 



(3.21) 



(3.22) 



These two ratios agrees well with the lattice calculations [0 which lead to 1.12 ± 0.05 and 
1.34 ± 0.07, respectively. 

The QCD corrections may be considered via the renormalization of the current J = qTQ 



II]. In renormalization-group-improved perturbation theory, up to the next-to-leading 



order the values of the decay and coupling constants at energy scale \x = tuq are given by 
fM(m Q ) 



A 



N T -F(m Q )[l + d M ^^-}[l + —(G 1 + 2d M G 2 )] 

\ m M A M 6tt m Q 



F{m Q ) 



, 1 + d M— R 1 + [9i + 2dM92)\, 



71 



2n 



(3.23) 
(3.24) 



with /5 = 33 — 2np, M = B, D, B*, D* and tuq = mi,, m c . From the results of /jy(2A) given 
in (|3.20p , we then have 

f B (m b ) = 0.159 ± 0.042 GeV, fB*{m b ) = 0.166 ± 0.038 GeV, 

f D (m c ) = 0.251 ± 0.099 GeV, f D , (m c ) = 0.293 ± 0.086 GeV. (3.25) 

which is consistent with the experimental upper limit < 200 MeV and /d(^d) < 

290MeV and also some theoretical upper bounds p| . The results also agree with the lattice 



calculations [|T3| and with full QCD calculations [18 
The ratio in eq.( |2.14j ) is now modified to be 

1/2 — 

f v m^ A P 1/2 2a s (m Q ) 



fpmj, 



1/2 



(1 



3tt 



_ 2a a (m Q ) 



3tt 



)(1 G 2 ) 

m Q 

> 

-#2), 



which yields 



/• 1/2 

fBm l B 2 
f 1/2 

jD* m D* 



1.05 ±0.03, 



1.22 ±0.07. 



(3.26) 

(3.27) 
(3.28) 
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IV. QCD SUM RULE EVALUATION ON WAVE FUNCTIONS f AND «j 

It is easily seen that to the leading order of heavy quark expansion, rriq — > oo, there is 
no difference between the HQEFT and the usual HQET, thus the procedure of calculating 
Isgur-Wise function £(?/) is the same as in the usual HQET. The main task in this section 
is to evaluate the two additional wave functions K\(y) and ^{y) involved at 1/rag order. 
For completeness, we also briefly outline the calculation of the Isgur-Wise functions £(y) in 
HQEFT. 

To evaluate the wave functions, we need to consider three point Green functions of the 
relevant operators. For the Isgur-Wise function £(y), it relates the following three-point 
correlation function at leading order of 1/rriQ 

E(u,u/) = J d\<£ye l{P '- x - p - y) < 0|(gf M ,Q\ x) , (QTQ) {0) , (QT M q) {y) \0 > 

= J tfixd^eW—**) < 0\(qf M/ Q'+) {xh (Q'+TQ+) {oh (Q+r M g) (2/) |0 > 

+ 0(l/m Q ), (4.1) 

where the Dirac structure T of the heavy-heavy current can in principle be arbitrary. For 
the present case, we take T = 7 M . As S(u;, to') is an analytic function in ui = 2v-k + 0(l/m,Q) 
and (J = 2v' ■ k' + 0(l/m,Q) with discontinuities on the positive real axis, one can write the 
phenomenological representation as 

0\qT M Q'+\M^ >< M' V ,\Q'+TQ+\M V >< M v \Q+r M q\0 > m M m M . 
A m A m ,(lj m -uj-ie) {u' M - u' - le) m z Q 

+ ( dvdp'- Pphya ^ V>) + subtractions. (4.2) 

J [y — to — it) [y 1 — uo 1 — it) 



The first term in (|4.2f) is a double-pole contribution and the second represents the higher 
resonance contributions in the form of a double dispersion integral over physical intermediate 
states. Parameterizing the three matrix elements in ( |4.2|) by using (|2.9|) and noticing the 
following relation 

Tr[t MI M(v')]Tr[M(v')TM(v)]Tr[M(v)T M ] = 4A 2 Tr[f M ,P' + rP + T M ], (4.3) 
the double-pole term becomes 



7 = Tr[T M P\TP+T M ] A 2 m M m M . z , 2 

pole (% — u — ie) (u' M — uj' — ie) Am^-m 1 rai 



F 2 £. (4.4) 



Theoretically, the correlation function (Q|) may be written as 



1 — 'theor 



dvdv' Ppert{ — ) _|_ w Np _|_ subtractions (4.5) 

[y — uj — it) \y' — lo — it) 



which can be calculated perturbatively in deep Euclidean region (u,u f A). Note that 
there are two momentum variables for the correlator (|4.1| ), a double Borel transformation 
^B&) should be applied to both sides of the sum rule. Because of the heavy quark 
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symmetry, 00 and u' are symmetric in ( f4.1|) , and thus it is natural and convenient to choose 
r = t' = 2T. It is useful to define u± = ± u'), and integrate first the spectral function 
over uo_ at the region — v + < v_ < u + . Finally the quark-hadron duality allows us to write 

E ph en = 2 dv + e u+/ p pert (v + ) +E NP ; (4.6) 

where H denotes the result obtained by applying double Borel operators to H, and 

Ppert{v+) = dv_p pert {v + ,V_). (4.7) 



The one loop perturbative diagrams and lowest order nonperturbative diagrams propor- 
tional to quark condensate and mixed quark-gluon condensate are listed in Fig. 2. In this 
section the gluon condensate can be safely neglected since its contribution is tiny. Calcula- 
tion of those Feynman diagrams in Fig. 2 gives 

-<qq>-i< iOvuF^q > (4.8) 

The continuum threshold energy in ( |4.8|) is in general a function of the recoil variable y. 
One may employ different models for reasonable choice of this function. It is seen that if 
the l/m,Q order terms and order a s terms in (|3.13|) are neg lected, flU) reduces to (gT3|) at 



the zero recoil point. This implies that we may use the same values of T and uj]\j as those 
in Q3.13 ) for evaluating the wave functions £, K\ and k 2 - The values of T and ujm can be 



read from Fig. 4, Fig. 5 and Fig.6. Note that the threshold energy satisfies the normalization 
co»o(l) = ujq. Our numerical results are plotted in Fig. 7, where we have used the values of 
ujq, and F obtained in the previous section. To be consistent, as the QCD radiative 
corrections are not considered in eq.( |4.8| ), we have used the values given in eq.( |3.18| ) where 
the results were obtained without QCD corrections. For comparison, we have used two 
simple models considered in 0: 

t \ /,\ f 1, model 1; . . 

^(»)=wb(l)|^ ( model2 . (4-9) 

We now turn to the calculations for the wave functions Ki(y) defined in ( |2.9| ). For that, 
one may consider the following three-point correlation function at I/tuq order: 



K(w,w') = J d 4 xd 4 ye i{p '- x - p - y) < 0\T{(qT M ,Q) { 



)(x), 

J d 4 xd 4 ye« k '- x - k ^ < 0\T{(qT M ,Q'+) {x) , 

^(0:tr-^_(^ ± ) 2 g+) ( o)> (Qtr M q) {y) }\o > +o(i/m 2 Q ) (4.io) 
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Saturating this three-point Green function with hadron states, one gets the double-pole 
contribution 



v~ l t\ m M m M 



pole uiqAm^m'^m -oj- ie)(uj' M - u' - ie) 
x < 0|gf M Q;t|M; >< M^\Q>+r^-?±-(tip ± ) 2 Qt\M v >< M v \Q+Tq\0 > . (4.11) 

ZTfln W ■ L) 



The heavy-heavy current in (|4.10|) contains both spin-symmetry conserving operator 



2^Q'^ T 1^d( D ^ 2 Qv and spin-symmetry violating operator ^Q'jT^ B P + ^a a pF a l 3 Q+. 
The hadronic matrix elements of the latter are parameterized by the wave functions n 2 (y) 
and ^3(2/). In this note, we may consider only the Feynman diagrams shown in Fig. 3, 
namely we will neglect radiative corrections as a first approximation. In such a treatment 
the resulting contributions from the spin-symmetry violating operator may proportional to 
the mixed quark-gluon condensate. Noticing that in the fixed point gauge 

< 0| : q(x)A li (z)q(0) : |0 > = Z —o V( , < qa a/3 F^q > 

+ higher dimensional condensates to be neglected, (4-12) 

it is readily seen from ( |2.12| ) that there would be no contributions to k 3 (?/) at the order we 
are considering. So in this approximation we have K^y) = 0, namely ^(y) only receives 
contributions from higher order and higher dimensional condensates which are expected to 
be small. For this reason, we may rewrite ( [4.1 1| ) as 

v 1 i\ Tr[T M ,P^TP + T M ] A 2 m M m M > „ 2 1 / , A \ / A ta\ 

lC pole {u, UJ ) = — r-^-j 2 F o + 4/K2 , (4.13) 

{uj M ~uj-ie){uj' M -uj' -ie) A M A M > m z Q 2m Q A 

where we have used (Of) and (J4.3[) as well as the formula P + a^ v A4(v )a^ u = 2cIm-M(.v). In 



evaluating the three-point Green function of ( 4.10|) , one meets a non-local operator, to be 



convenient of calculating /C, one may choose the axial gauge v ■ A — 0. Note that in this 
gauge the diagrams with gluons attached to a heavy quark line are absent. 

Adopting the same strategy for evaluating the Isgur-Wise function £(y), we arrive at the 
following double Borel transformed sum rule 

2m Q A 2n 2 Jo m Q{y + 1) 

By separately considering the spin-conserving and spin-breaking corrections to the limit 
case mq — ► 00 (or equivalently considering those terms with and without (1m), and taking the 
central values in fl3.18|) , we obtain numerical results which are plotted in Fig.8. The curves 
in Fig.8 are corresponding to the results at T = l.OGeV which is chosen by considering the 
stability region of the curves exhibited in Fig.4. We also present in Fig. 9 amd Fig. 10 the 
values of ^i(l) and k 2 (1) as functions of the Borel parameter T. It is seen from those two 
figures that ^i(l) and ^(1) are really stable at the region around T = l.OGeV. The stable 
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regions in Fig. 9 and Fig. 10 are consistent each other and also consistent with the ones 
in Fig.4-Fig.6. With these considerations and analyses, we may present the final result for 
«i(l) as 

Kx = Kl (l) = -0.50±0.18GeV 2 (4.15) 

which agrees with the one extracted from the heavy meson masses 0, where «i(l) could 
range from — 0.8GeV 2 to — 0.25GeV 2 with a favoriable value «i(l) ~ —0.6lGeV 2 . 

In our previous papers fH-§J, we have argued that < iv ■ D > is of order the binding 
energy A. For simplifying the analyses, we have actually made a heavy quark expansion at 
point < iv ■ D >= A for inclusive decays of heavy hadrons To check the validity of this 
approximation, we may replace the non-local operator in ( |2.9| ) with 4 and evaluate the 



resulting local matrix element, then a comparison between two results should allow one to 
test the goodness of the approximation. By doing this, we may reparameterize the matrix 
elements as 

< M^\Q'+T^P + DlQ+\M v >= -Ki(y)^Tr[M'TM], 

< M^\Q'+T^P + l -a a ^Q+\M v >= lTr[K a ^v,v')MTP + l -a^M], (4.16) 

where K a p may be decomposed in a similar way as K a p. Applying the sum rule approach 
once more to evaluate the parameter K\ and following the same strategy as before, we yield 
the following sum rule formula for K\, 

where the corresponding Feynman diagrams are the same as those in Fig. 3 with the box 



now representing the new operator — ^^Q'^T~-(iIp^) 2 Q^ . The numerical results of ([4. 17]) 
are shown in Fig. 11 and Fig. 12, it is found that 

Kx(l) w -0.40GeV 2 (4.18) 

which is slightly lower than K\(l). Comparing the two nemerical results of K\(l) and «i(l), 
we see from ( |2.9| ) and (|4.16|) that the simple replacement 

iv-D~vkr*iA (4.19) 

is actually a reliable approximation for the operator iv ■ D. 

The above demonstration supports the analyses in We are confirmed to believe 

that the HQEFT is more reliable to describe the off-mass shell heavy quark within heavy 
hadrons and can provide a consistent understanding on both exclusive and inclusion decays 
of heavy hadrons. 

We may return to comment on the wave function K2(y). Up to the order considered 
above, the sum rule formula Q4.14 ) leads to the value ^2(1) ~ 0.015GeV 2 which is much 



smaller than the one extracted from the meson masses [0 (where k 2 (1) w 0.056GeV 2 ). One 
of possible reasons for such a discrepancy may be seen from ( |4.14|) in which only one mixed 
condensate term contributes to k 2 . That unique term is relevant to the spin-symmetry 
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breaking term and arises from the diagram Fig.3(c). As the operator parameterized by 
k 2 in Q2.9Q contains gluon fields, the lowest order perturbative contributions should arise 
from the two- loop diagram in Fig. 11. In order to improve the determination for k 2 , one 
should therefore calculate at least two-loop perturbative diagrams as well as nonperturbative 
diagrams up to order a s in a consistent way. Since all these diagrams are at least at the 
order of a s , the more precise value of n 2 should not be too large. We should not present 
calculations of the additional diagrams for k 2 in this paper. 

In fact there is another way to estimate the values of «i(l) and n 2 (l) directly from the 
results obtained through two-point correlator functions in Sec.[TTI|. This is again because the 
remarkable relation between the heavy-light hadron mass and wave functions at zero recoil 
resulted from the complete HQEFT. It can be seen from Eqs. (|2.15|) and ( |3.14j ) that the form 
factors fti(l) and ^2(1) are simply given by 

«i(l) = -^f-, - 2 (1) = -^- (4.20) 



With these relations, we obtain from the values of lu^ in eq. (|3.19|) 



«i(l) « -0.43GeV 2 , k 2 (1) « 0.08GeV 2 , (4.21) 



which are consistent with the results yielded above (eq. (|4.15 )) and also those obtained in 
Ref. M for «i(l) and k 2 (1). 



V. COORECTIONS FROM TWO-LOOP PERTURBATIVE QCD 

In order to take a look at the magnitude of the effects of QCD radiative corrections, as 
in Ref. ||, one can now include the two- loop perturbative contributions. Their effects can 
be simply taken into account by replacing the perturbative contributions in the sum rule 
( P. 11 ) with the following ones 



tt t s 3 r c , t 1 1 2/1 2a, ri 2A 13 27r 2 n 3u , , ri , 

IWM = — / dv{- + -)v\\ + —Z[ln— + — + — - ). (5.1) 

6ir z Jo 4mQ v — uo — le n v o 9 4mg 

In comparison with the results from leading QCD corrections, we have plotted the modified 
results with two-loop QCD corrections in Figs. (14-22). As a consequence, instead of 
eqs. fl3.18| ) and (|3.19| ), we arrive at the following modified results 



and 



uj = 1.8 ± 0.3 GeV, 

A = 0.56 ± 0.08 GeV, 
2 

F = 0.38 ± 0.06 GeV 3/2 . (5.2) 



uj x = 1.0 ± 0.2 GeV 2 , 
cjW = 0.65 ±0.10 GeV 2 , 



16 



G x = 


0.75 ±0.15 GeV, 


9i = 


0.46 ±0.12 GeV, 


LU 2 = 


-0.15 ±0.05 GeV 2 , 


, ; (2) - 
UJ M — 


-0.14 ±0.03 GeV 2 , 


G 2 = 


-0.09 ±0.03 GeV, 


92 = 


-0.06 ±0.02 GeV. 



(5.3) 



Correspondingly, instead of eq. (|3.25 ), we have 



f B (m b ) = 0.196 ± 0.044 GeV, 
f D (m c ) = 0.298 ± 0.109 GeV, 



f B *{m b ) = 0.206 ± 0.039 GeV, 
f D *{m c ) = 0.354 ± 0.090 GeV. 



(5.4) 



Comparing these values with those obtained in Sec.|ITT|, we see that the QCD radiative 
corrections may enlarge F by about 25%. It implies that the radiative corrections may be 
significant for a more accurate determination for some physical quantities. The large values 
seems to be consistent with the ones from the recent calculations by Lattice QCD approach 



IMZI- 

With those values in eq.i 



2|), we yield from the sum rules in eqs. (|4.14|) and ( [4.1 7|) 



-0.34GeV 2 ; 
-0.26GeV 2 . 



(5.5) 
(5.6) 



which are lower than the ones without two-loop perturbative QCD corrections. While we 
would like to point out that the modified results for /«i(l) and Ki{l) ( eqs.(|5.5[) and ( |5.6| ) ) 
may not be regarded to be more reliable than the ones given in eqs.( [4.15"|) and ( [4.1 8|) since 
eqs.( 4.14| ) and ( 4.17 ) contain only the leading order contributions in perturbation theory. 
For a complete and consistent evaluation, one should also include the next-to-leading order 
corrections to eqs.( |4.14| ) and ( [4.17| ). But comparing the two values in eqs.( |5.5| ) and ( |5.6| ) we 
see that the approxmation iv ■ D ~ v ■ k ~ A still holds, though both the values of ^i(l) 
and Ki(l) are now smaller than those presented in Sec.fVJ because the input parameter F 
is enlarged by the two-loop perturbative corrections. 

It is also interesting to notice that the relation Ki(l) = ~ — 0.37GeV 2 is still 



satisfied well by looking at the result given in eq.(|5.5|) and the value of in eq 
conclusion, all the numerical results in this paper turn out to be consistent. 



In 



VI. CONCLUSIONS AND REMARKS 



We have present, within the complete HQEFT, a consistent evaluation for the decay 
constants and wave functions of heavy- light hadron systems up to order of l/mq. It has 
been seen that the QCD Lagrangian for heavy quarks do neet to be transformed into a new 
heavy quark effective Lagrangian with including contributions of both particle fields and 
antiparticle fields 0, and the currents containing heavy quark can be consistently expanded 
in powers of l/mq. Though the leading operator in the l/mq expansion is (must be) the 
same as the one in the usual HQET, the operators at order 1/mg begin to be different from 
those in the usual HQET. 
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In the complete HQEFT, corrections arising from both the current expansion and the 
insertion of Lagrangian into the heavy-light meson anahilation matrix elements are charac- 
terized by two factors G\ and G 2 . Similarly, corrections arising from both sources to the 
heavy meson transition matrix elements are characterized by three functions K\(y), n 2 (y) 
and Ks(y). Furthermore, the values of k,\ and k 2 at the zero recoil point also characterize 
the 1/rriQ order corrections to the meson masses. These makes the HQEFT be much elegant 
than the usual HQET. 

The leading order contributions and the 1/mq corrections to heavy meson decay con- 
stants and heavy meson transition matrix elements have been investigated consistently by 
using QCD sum rule approach within the framework of HQEFT of QCD. For heavy-light 
meson decays, we have calculated the form factor F at leading order (tuq — > oo), and the 
form factors G\ and G 2 concerned at the 1/tuq order as well as the binding energy A. Par- 
ticularly, we have found that the 1/mQ corrections to the heavy-light meson decay constants 
are actually determined by two composite form factors gi and g 2 . These two composited 
form factors were found to be much smaller than the heavy quark masses, which implies 
that the scaling law of the decay constants are only sligtly breakdown. This observation 
shows that the 1/rriQ expansion in the conplete HQEFT works well, which is unlike the 
usual HQET which may lead to, as shown in [Q, the breakdown of the 1/rriQ expansion in 
evaluating the decay constants. Our results for the heavy-light meson decay constants have 
also shown a good agreement with the known experimental results and upper limits. 

We have also calculated the Isgur-Wise function and 1/rriQ order spin-symmetry con- 
serving form factor K\ as functions of the recoil value. It have been found that ~ 
-0.50±0.18GeV 2 , which agrees with the value extracted from the interesting relations be- 
tween meson mass and wave functions at zero recoil. We have also illustrated how the simple 
replacement iv ■ D ~ v ■ k rs A holds. This shows that the residual momentum of heavy 
quark within heavy-light hadron does be around the binding energy which had been seen 
H to be the main point to understand the puzzle of the bottom hadron lifetime differences. 

The spin symmetry breaking factor k 2 has been discussed. The only diagram at leading 
order yields a k 2 value which is much smaller than that obtained in ||. Thought its value has 
not yet been evaluated accuratly, it implies that k 2 must be small as it characterizes the spin 
symmetry breaking effects of heavy-light mesons. A further calculation of k 2 including higher 
order contributions remains an interesting subject. ^{y) was also found to be small since 
it receives contributions only from higher order radiative diagrams and higher dimensional 
condensates. 

An interesting feature resulting from the HQEFT is that the values of fti(l) and k 2 {1) 
can also be simply obtained, due to the interesting relation between heavy-light meson mass 
and wave functions at zero recoil, from two-point Green's function in evaluating the heavy- 
light meson decay constants via sum rule approach. It is remarkable that the resulting 
values of ^i(l) and k 2 {1) in this way do agree with the results obtained from the analysis 
of the three-point Green's function. Finally, we have shown that the higher order radiative 
corrections could be nontrival for a more accurate calculation of heavy-light meson decays. 
But the relations iv ■ D ~ iv ■ k ~ A and fti(l) = — a% /2, k 2 (1) = — a% /2 hold even when 
higher radiative corrections are included. 

In summary, the complete HQEFT works well for describing the slightly off-mass shell 
heavy quark within heavy hadrons. In this paper, we have further checked the consistent 
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of the HQEFT in applying to the heavy-light hadron systems with tuq » A » m q . For 
heavy-heavy bound state systems, such as bottom-charm system like B c , and muonium (/j,e) 
system in QED, one needs to make 1/M expansion for both heavy particles and to redefine 
the effective fields and bound states when applying for the complete heavy particle effective 
field theory with keeping the antiparticle contributions ||TJ. 
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FIGURES 




Fig.l. Feynman diagrams contributing to heavy meson decay. The thick lines are heavy 
quarks; the light lines are light quarks; the curves are gluon fields; the black dots represent 
condensates; and the external dashed lines are the heavy-light currents considered in ( |3.1| ). 
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Fig. 2. The lowest order Feynman diagrams contributing to £. The wave lines represent the 
heavy- heavy current Q'^TQ^ in (p~l|). 








Fig. 3. The lowest order Feynman diagrams contributing to k% and «2- The box at the up of 
each diagram represent the 1/mg order heavy- heavy current — 2Wi~Q'&^l^T) (^-i-) 2 Qt m ( 4.10Q . 
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Fig.4. Leading oreder sum rule result for heavy meson decay. 
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Fig.5c 

.5 Sum rule result for 1/mg order spin-symmetry conserving corrections to heavy meson decay. 
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Fig.6b 
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Fig.6c 

Fig.6. Sum rule result for I/tuq order spin-symmetry breaking corrections to heavy meson decay. 
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Fig. 7. Sum rule result for Isgur-Wise function £,(y). 
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model 1 
model 2 
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Fig.8. Sum rule result for l/mq order transition form factor K\ and ^2- 
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Fig.9. ^i(l) as a function of Borel parameter T. 
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Fig.10. «2(1) as a function of Borel parameter T. 
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Fig. 11. Sum rule result for K\. 
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Fig.12. K\{\) as a function of Borel parameter T. 




Fig. 13. The lowest order perturbative diagram contributing to «2. 
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Fig. 14. Leading oreder sum rule result for heavy meson decay when two- loop perturbative contri- 
butions are considered. 
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Fig. 5 Sum rule result for 1/m.Q order spin-symmetry conserving corrections to heavy meson decay 
when two-loop perturbative contributions are considered. 
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Fig. 16c 

Fig. 6. Sum rule result for 1/tuq order spin-symmetry breaking corrections to heavy meson decay 
when two-loop perturbative contributions are considered. 




33 



nodel 1 
model 2 



.6 1 .7 1 .8 1 .9 



Fig. 18. Sum rule result for 1/uiq order transition form factor K\ and K2 when the value of F in 
eq.([5.2D is used as an input parameter. 
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Fig. 19. K]_(l) as a function of Borel parameter T when the value of F in eq.(p 
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Fig. 20. ^2(1) as a function of Borel parameter T when the value of F in eq.( |5.2| ) is used as an 
input parameter. 



model 1 
model 2 



.2 1 .3 



1 .5 1.6 1 .7 



Fig.21. Sum rule result for K\ when the value of F in eq.(5.2) is used as an input parameter 
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Fig. 22. ^i(l) as a function of Borel parameter T when the value of F in eq. fl5.2D is used as an 
input parameter. 
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